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Abstract 

A Green's function approach to the inclusive quasielastic (e, e') scattering 
is presented. The components of the nuclear response are written in terms 
of the single-particle optical model Green's function. The explicit calculation 
of the Green's function can be avoided by its spectral representation, which 
is based on a biorthogonal expansion in terms of the eigcnfunctions of the 
non-Hcrmitian optical potential and of its Hermitian conjugate. This allows 
one to treat final state interactions consistently in the inclusive (e, e') and in 
the exclusive (e,e'N) reactions. Numerical results for the longitudinal and 
transverse response functions obtained in a nonrelativistic and in a relativistic 
framework are presented and discussed also in comparison with data. 



1 Introduction 

Inelastic electron scattering has proven to be a very precise tool for studying the 
excitation spectrum of atomic nuclei owing to the possibility of simultaneously 
varying energy and momentum transfer (w, q) [1]. In a representation of the nuclear 
response as a function of ui and q a large broad peak occurs at about u = q 2 /2M. Its 
position corresponds to the elastic peak in electron scattering by a free nucleon. It is 
quite natural to assume that a quasifree process is reponsible for such a peak with a 
nucleon emitted quasielastically. Coincidence (e, e'N) experiments in the quasifree 
region confirm such a picture and represent a valuable source of information on 
single- nucleon degrees of freedom inside nuclei. 

In the inclusive (e, e') process only the scattered electron is detected and the 
final nuclear state is undetermined, but the main contribution in the region of 
the quasielastic (QE) peak comes from the interaction on single nucleons. If the 
nucleons were indeed free, the peak would be sharp and would just occur at the 
energy loss u = q 2 /2M corresponding to the energy taken by the recoiling free 
nucleon. A shift in the position of the peak is produced by the nuclear binding, 
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while a broadening of the peak is produced by Fermi motion. A very simple model, 
where the cross section is given by the interaction on all the constituent nucleons of 
the nucleus and a Fermi gas model is assumed for nuclear structure, is indeed able 
to give a remarkably good description of the cross section in the QE region [2]. This 
same model, however, is not able to describe simultaneously the QE longitudinal 
and transverse responses, whose Rosenbluth separation has been experimentally 
achieved on a variety of nuclei and over a large range of momentum transfers [1]. 

In the one-photon exchange approximation and assuming the plane-wave 
approximation for the incident and the outgoing electrons, the (e, e') cross section 
is given by [1] 

a inc = K (2e L R L + R T ) , (1) 
where K is a kinematical factor and 

£L = ^l + 2^tan 2 (V2)j (2) 

measures the polarization of the virtual photon. In Eq. (2) i? e is the scattering 
angle of the electron, Q 2 = q 2 —uj 2 , and q^ = (to, q) is the four momentum transfer. 
The longitudinal and transverse response functions, R\, and Rt are defined by 

R L (u>, q) = W&(<J, q), Rr(u, q) = W%(u, q) + W&(u, q) , (3) 

in terms of the diagonal components of the hadron tensor 

WZip, Q) = I £ I <*f I JM (9) I *o> I 2 S(E +u-Ef). (4) 
f 

Here J M is the nuclear charge-current operator which connects the initial state 
| ^o) of the nucleus, of energy Eq, with the final states | of energy Ef. In the 
inclusive (e, e') process the sum is over a complete set of final nuclear states and 
only the diagonal components of the hadron tensor contribute, while in the exclusive 
(e, e'N) process also the nondiagonal components contribute and the cross section 
is in general given in terms of a larger number of structure functions [1]. 

The separation of the response functions R^ and i?x can experimentally by 
achieved from Eq. (1) by varying electron kinematics. In a plot of the cross section 
as function e^, for fixed values of uj and Q 2 , the slope gives R^ and the extrapolated 
intercept with the vertical axis at £l = gives Rt. Such a Rosenbluth separation 
[3] has been achieved in the 80's and 90's on different nuclei [1]. 

A huge amount of theoretical work was produced over the past two decades 
in order to explain the problem raised by the experimental separation of R^ and 
i?T [1]- The main problem was the apparent quenching of i?L, while for R^ there 
is in general an apparent excess of strength. Different approaches, with different 
theoretical ingredients, obtained partial success in explaining either the longitudinal 
or the transverse response, but despite all these efforts some problems are still 
open and a consistent and simultaneous description of i?L and Rt data has not 
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been achieved. Some confusion is also due to the experimental situation. These 
experiments of separation are difficult and data from different laboratories are 
sometimes in disagreement. New experiments with improved experimental accuracy 
would be helpful to make the experimental situation clearer. New measurements 
are planned at Jlab [4] 

In this contribution we present a Green's function approach to the inclusive 
QE (e, e') scattering. With our work we don't aim at a unified and consistent 
description of i?L and Rt data. The model contains approximations. It is basically 
a single-particle approach, where only the one-body part of the nuclear current 
is retained. Calculations of two-body contributions performed by different groups 
have given different results [1, 5]. There are, however, consistent indications that the 
combined effect of two-body current and tensor correlations can be appreciable on 
the transverse response [6, 7, 5]. Thus, with our model we don't expect to describe 
Rt data, while the main contributions to Rl should already be included. Our aim is 
to study the effects of final state interactions (FSI). FSI are an important ingredient 
of the inclusive electron scattering, since they are essential to explain the exclusive 
one-nucleon emission, which gives the main contribution to the inclusive reaction in 
the QE region. The absorption in a given final state due, e.g., to the imaginary part 
of the optical potential, produces a loss of flux that is appropriate for the exclusive 
process, but inconsistent for the inclusive one, where all the allowed final channels 
contribute and the total flux must be conserved. 

This conservation is preserved in the Green's function approach, where the 
components of the nuclear response are written in terms of the single-particle optical 
model Green's function. This result was originally derived by arguments based 
on the multiple scattering theory [8] and successively by means of the Feshbach 
projection operator formalism [9, 10, II, 12]. The spectral representation of the 
single-particle Green's function, based on a biorthogonal expansion in terms of 
the eigenfunctions of the non-Hermitian optical potential, allows one to perform 
explicit calculations and to treat FSI consistently in the inclusive and in the exclusive 
reactions. Important and peculiar effects are given, in the inclusive reactions, by the 
imaginary part of the optical potential, which is responsible for the redistribution 
of the strength among different channels. The approach has been developed and 
used in a nonrelativistic frame [11] and, more recently, also in a relativistic frame 
[13]. Although some differences and complications are due to the Dirac matrix 
structure, the formalism follows in both frames the same steps and approximations. 
The numerical results obtained in both frames allow us to check the relevance of 
relativistic effects. 

The Green's function approach is presented in Sec. 2. In Sec. 3 some results for 
i?L and Rt are presented and discussed also in comparison with data. Summary 
and conclusions are drawn in Sec. 4. 
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2 The Green's function approach 



The response functions of the inclusive (e, e') scattering are defined in Eq. (3) 
in terms of the diagonal components of the hadron tensor of Eq. (4), that can 
equivalently be expressed as 

W&(u,q) = --Im(* I jrt(q)G(Et)Jl>(q) | * ) , (5) 

where Ef = Eq + uj and G{E{) is the Green's function related to the nuclear 
Hamiltonian H, i.e., 

°< £ <> " e^hTV, ' < 6 > 

Here and in all the equations involving G the limit for 7] — > +0 is understood. 
It must be performed after calculating the matrix elements between normalizable 
states. 

The hadronic tensor in Eq. (5) contains the full (A+l)-body propagator of the 
nuclear system and, as such, it is an extremely complicated object, which defies a 
practical evaluation. Only an approximated treatment reduces the problem to a 
tractable form. Under suitable approximations the nuclear response can be written 
in terms of the optical-model Green's function [11]. 

The first approximation consists in retaining only the one-body part of the 
charge-current operator J' 4 . Thus, we set 

A+l 

^(<7)=£jf(<7), (7) 

i=i 

where jf acts only on the variables of the nucleon i. By Eq.(7), one can express the 
hadron tensor as the sum of two terms, i.e., 

W&(u,, q) = W*( U , q) + W&(u,, q) , (8) 

where W^^(uj,q) is the incoherent hadron tensor [14], which contains only the 
diagonal contributions jf^Gjf, whereas the coherent hadron tensor W^(tu,q) 
gathers the residual terms of interference between different nucleons. As the 
incoherent hadron tensor, also W^(u, q) can be expressed in terms of single-particle 
quantities (see Sect. 9 of Ref. [12]), but for the transferred momenta considered here 
we can take advantage of the high-g approximation [15] and retain only W^(uj, q). 
This term can be further simplified using the symmetry of G for the exchange of 
nucleons and the antisymmetrization of | \£o)- Therefore, we write 

W&(u,q) ~ W^(uj,q) = -^-±Alm<* | j^(q)G(E { )f (q) \ * ) , (9) 

7T 

where j^{q) is the component of J^{q) related to an arbitrarily selected nucleon. 

Let | n) and | e) denote the eigenvectors of the residual Hamiltonian Hr of 
A interacting nucleons related to the discrete and continuous eigenvalues e n and e, 
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respectively. We introduce the operators P n , projecting onto the n-channel subspace 
of 7i, and Q n , projecting onto the orthogonal complementary subspace, i.e., 

P n = [ dr\ r;n){n;r |, Q n = l-P n . (10) 

Here | r;n) is the vector obtained from the tensor product between the discrete 
eigenstate | n) of Hr, and the orthonormalized eigenvectors | r) of the selected 
nucleon. Moreover, we introduce the projection operator onto the continuous 
channel subspace, i.e., 

P c = J de J dr | r;e){e;r | . (11) 

Due to the completeness of the set {| r;n),\ r;e)}, one has 

£P n + P c = l. (12) 

n 

Then, we insert Eq. (12) into Eq. (9) disregarding the contribution of P c . This 
approximation, which simplifies the calculations, is correct for sufficiently high 
values of the transferred momentum q. Thus, the hadron tensor of Eq. (9) can 
be expressed as the sum 

WW(u, q) = O, q) + Wg£(u, q) , (13) 

of a direct term 

<>,<?) = £wr (",?), 

n 

WW{u,q) = -^±llm(* | f\q)P n G{E i )P n f{q) \ * ) , (14) 

7T 

and of a term 

WZfaq) = ^Wr(co,q), 

n 

Wr(oJ,q) = -^±±hn(* | f\q)P n G{E,)Q n f{q) \ * ) , (15) 

which gathers the contributions due to the interference between the intermediate 
states | r;n) related to different channels. 

If we disregard the effects of interference between different channels and consider 
only the direct contribution to the hadron tensor of Eq. (14), the matrix elements 
of P n G(E)P n in the basis | r;n) define the single-particle Green's function Q n (E) 
[11] of the single-particle Feshbach optical potential V n (E) [17] 

e " (£) " E -T-V„( E) + i„ ■ (16> 
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which describes the elastic scattering of a nucleoli by an A-nucleus in the discrete 
state | n). Using the definition of P n in Eq. (10), the direct hadron tensor W^(uj, q) 
of Eq. (14) can be reduced to the single-particle expression 

W{?(u>, q) = --A n Im(^ n | jrt(q)G n (Et - e n )f(q) \ <p n ) , (17) 

7T 

where A n is the spectral strength [16] of the hole state <p n , which is the normalized 
overlap integral between | *$>o) and | n). The problem of expressing the interference 
hadron tensor in a one-body form is treated in Ref. [11]. It is argued that 

the contribution of can be included into the direct hadron tensor by the 
simple replacement 

g n (E) - gf(E) ee ^l-V n (E)G n {E)^l-V n (E) , (18) 

where V n {E) is the energy derivative of the Feshbach optical potential. 

In principle, different optical potentials V n must be considered for different 
values of n. As neither microscopic nor empirical calculations are available for the 
optical potential V n associated with the excited states, a common practice relates 
them to the ground state potential Vo by means of an appropriate energy shift. 
Therefore we set 

V n (E) ~ V (E) , (19) 

which implies 

g n (E)~g (E). (20) 

Using these approximations, we write 

WW(u,q) = --^A n Inu>„ | fHq)gf(E i - e n )f{q) \ <p n ) . (21) 

As a next step, the spectral representation of the single-particle Green's function 
function can be used in order to allow practical calculations of the hadron tensor of 
Eq. (21). 

In expanded form, the hadron tensor reads 

W»(u,q)= - -^2^nlm( l p n \j^(q)Jl-V(E) 

x g(E)^jl-V{E)f(q) | <p n ) , (22) 

where E = Ef — e n . Here and below, the lower index is omitted in the Green's 
functions and in the related quantities. 

Due to the complex nature of V(E) the spectral representation of Q(E) involves 
a biorthogonal expansion in terms of the eigenfunctions of V(E) and of its Hermitian 
conjugate V\E). We consider the incoming wave scattering solutions of the 
eigenvalue equations, i.e., 

(S-T-V\E))\ X ^(E)) = 0, (23) 

(£-T-V(E))\x { £ -\E)) = 0. (24) 
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The choice of incoming wave solutions is not strictly necessary, but it is convenient 
in order to have a closer comparison with the treatment of the exclusive (e, e'N) 
reactions, where the final states fulfill this asymptotic condition and are the 
eigenfunctions | \~e\e)) of V\E). 

The eigenfunctions of Eqs. (23) and (24) satisfy the biorthogonality condition 



(x { i\E)\x { £ -\E)) = 5 (£-£') , 
and, in absence of bound eigenstates, the completeness relation 

H &£\x { s\e))(x [ s\e)\=1 
Jm 



(25) 



(26) 



where the nucleon mass M is the threshold of the continuum of the Feshbach 
Hamiltonian. 

Eqs. (25) and (26) are the mathematical basis for the biorthogonal expansions. 
The contribution of possible bound state solutions of Eqs. (23) and (24) can be 
disregarded in Eq. (26) since their effect on the hadron tensor is negligible at the 
energy and momentum transfers considered in this paper. 

Using Eqs. (26) and (24), one obtains the spectral representation 



Therefore, Eq. (22) can be written as 

1 r r°° 

WW(u,q) = Vim / d£ 



Et 



£ + 177 



TW{£,Et-e n ) 



(27) 



(28) 



where 



Tr(£,E) = X n ^ n \f\q)^l-V{E)\x { £ \E)) 
x {x£~\e) \\Jl- V{E)f(q) I ip n ) . 



(29) 



The limit for r] — > +0, understood before the integral of Eq. (28), can be calculated 
exploiting the standard symbolic relation 



lim 



1 



V 



1 



ivrJ (E - £) 



n-^oE-S + 'v} \E-£ 
where V denotes the principal value of the integral. Therefore, Eq. (28) reads 



(30) 



W^(uJ,q) 



E 

n 

w 

7T JA 



ReT^{Ef-£ n ,Et-e n 
1 



d£- 

m Ef — e n — £ 



hnT^{6,Et-e n ) 



(31) 
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which separately involves the real and imaginary parts of . 

Some remarks on Eqs. (29) and (31) are in order. Disregarding the square 
root correction due to interference effects, one observes that in Eq. (29) the second 
matrix element (with the inclusion of \f\^) is the transition amplitude for the single 
nucleon knockout from a nucleus in the state | ^o) leaving the residual nucleus in 
the state | n). The attenuation of its strength, mathematically due to the imaginary 
part of V^, is related to the flux lost towards the channels different from n. In the 
inclusive response this attenuation must be compensated by a corresponding gain 
due to the flux lost, towards the channel n, by the other final states asymptotically 
originated by the channels different from n. In the description provided by the 
spectral representation of Eq. (31), the compensation is performed by the first 
matrix element in the right hand side of Eq. (29), where the imaginary part of V 
has the effect of increasing the strength. Similar considerations can be made, on 
the purely mathematical ground, for the integral of Eq. (31), where the amplitudes 
involved in have no evident physical meaning as £ / Ef — e n . We want to stress 
here that in the Green's function approach it is just the imaginary part of V which 
accounts for the redistribution of the strength among different channels. 

The matrix elements in Eq. (29) contain the mean field V(E) and its Hermitian 
conjugate V\E), which are nonlocal operator with a possibly complicated matrix 
structure. Neither microscopic nor empirical calculations of V(E) are available. 
In contrast, phenomenological optical potentials are available. They are obtained 
from fits to experimental data, are local and involve scalar and vector components 
only. The necessary replacement of the mean field by the empirical optical model 
potential is, however, a delicate step. 

In the nonrelativistic treatment of Refs. [11, 18] this replacement is justified on 
the basis of the approximated equation (holding for every state | ip) ) 

Im(V> | yfl-V(E)g(E)yfl-V'(E) | V) 
~ Im(V | yJl-V L (E)g L (E)yJl-V' L (E) \ , (32) 

where Vl(E) is the local phase-equivalent potential identified with the phenomeno- 
logical optical model potential and Gl(E) is the related Green's function. We have 
reasonable confidence that Eq. (32) holds also in the relativistic context. 

3 Results 

The response functions of the inclusive QE electron scattering are calculated from 
the single-particle expression of the coherent hadron tensor in Eq. (31). After 
the replacement of the mean field V(E) by the empirical optical model potential 
Vl(E), the matrix elements of the nuclear current operator in Eq. (29) are of 
the same kind as those giving the transition amplitudes of the electron induced 
nucleon knockout reaction in the distorted wave impulse approximation (DWIA) 
[1]. Thus, calculations have been performed adopting the same relativistic [20] and 
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nonrelativistic treatments [19] which were successfully applied to describe exclusive 
(e, e'p) data, and with the same phenomenological ingredients for bound state wave 
functions and optical potentials. 

In the calculations the residual nucleus states | n) are restricted to be single 
particle one-hole states in the target. A pure shell model is assumed for the nuclear 
structure, i.e., we take a unitary spectral strength for each single particle state and 
the sum runs over all the occupied states. In this way we are able to account for 
the contributions of all the nucleons in the nucleus, but correlations are completely 
neglected. This is of course an approximation that anyhow allows us to perform 
relatively simple calculations on a conceptuallly clear basis. 

The hadron tensor in Eq. (31) is the sum of two terms. The calculation of the 
second term, which requires integration over all the eigenfunctions of the continuum 
spectrum of the optical potential, is a rather complicate task. The contribution of 
this term is very small and can be neglected in the nonrelativistic frame [11], while 
in the relativistic frame it can be significant and must included in the calculation 
[13]. 

The longitudinal and transverse response functions calculated in the relativistic 
frame for 12 C at q = 400 MeV/c are displayed in Fig. 1 and compared with the 
Saclay data [21]. The low energy transfer values are not given because the relativistic 
optical potentials are not available at low energies. 

The agreement with the data is generally satisfactory for the i?L- In contrast, Rt 
is underestimated. This is a systematic result of both relativistic and nonrelativistic 
calculations. It may be attributed to physical effects which are not considered 
in the present approach, e.g., two-body currents. The effect of the integral in 
Eq. (31) is also displayed. At variance with the nonrelativistic result [11], this 
contribution is important and essential to reproduce the experimental longitudinal 
response. The contribution of interference between different channels, that produces 

the factor \Jl — V'^(E), gives only a negligible contribution in Fig. 1, while in the 
nonrelativistic calculation this contribution is important to improve the agreement 
with i?L data [11]. 

The contribution from all the integrated single nucleon knockout channels is 
also drawn in Fig. 1. It is significantly smaller than the complete calculation. The 
reduction, which is larger at lower values of to, gives an indication of the relevance 
of inelastic channels. 

An example of the comparison between the results of the relativistic and the 
nonrelativistic approaches is presented in Fig. 2 for the longitudinal response at 
q = 400 MeV/c. Both complete calculations are in satisfactory agreement with the 
data. This result is however due to a different effect of the v arious cont ributions in 
the two approaches. In the nonrelativistic case the factor \Jl — V L (E) produces a 
substantial reduction of the calculated response, that is necessary to reproduce the 
data. The integral in the second term of Eq. (31) gives only a small contribution and 
is neglected in the calculation. In the relativistic approach the integral is essential 
to reproduce the data while the interference between different channels is generally 
negligible. 
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The longitudinal and transverse response functions calculated in the relativistic 
approach for 12 C at q = 500 are displayed in Fig. 3 and compared with the Saclay 
data [21]. As already found in Fig. 1 at q = 400 MeV/c, a good agreement with the 
data is obtained for i?L, while i?x is underestimated. Also in Fig. 3 only a slight 
effect is given by the interference between different channels. The role of the integral 
in Eq. (31) decreases increasing the momentum transfer. The effect of the inelastic 
channels, indicated in the figure by the difference between the complete results and 
the contribution from all the integrated single nucleon knockout channels, is always 
visible and even sizable, but it decreases increasing the momentum transfer. 

The response functions for 40 Ca at q = 450 MeV/c are shown in Fig. 4 and 
compared with the available data. The calculated results are of the same order of 
magnitude as the MIT-Bates data, while for the Saclay data i?L is overestimated 
and i?T underestimated. The factor \Jl — V[JE) produces and enhancement which 
is minimal but visible in the figure. 

4 Conclusions 

A Green's function approach to inclusive QE electron scattering has been presented. 
The components of the hadron tensor are written in terms of Green's functions of the 
optical potentials related to the various reaction channels. The projection operator 
formalism is used to derive this result. An explicit calculation of the single-particle 
Green's function can be avoided by means of its spectral representation, based 
on a biorthogonal expansion in terms of the eigenfunctions of the non-Hermitian 
optical potential V(E) and of its Hermitian conjugate. The interference between 
different channels is taken into account by the factor yT — V'{E), which also allows 
the replacement of the mean field V{E) by the phenomenological optical potential 
Vl(-E'). After this replacement, the nuclear response functions are expressed in 
terms of matrix elements similar to the ones which appear in the exclusive one 
nucleon knockout reactions, and the same DWIA treatment can be applied to the 
calculation of the inclusive electron scattering. 

The effects of FSI are thus described consistently in exclusive and inclusive 
processes. Both the real and imaginary parts of the optical potential must be 
included. In the exclusive reaction the imaginary part accounts for the flux lost 
towards other final states. In the inclusive reaction, where all the final states are 
included, the imaginary part accounts for the redistribution of the strength among 
the different channels. 

All the final states contributing to the inclusive reaction are contained in the 
Green's function, and not only one nucleon emission. Our calculations for the 
inclusive electron scattering are different from the contribution of integrated single 
nucleon knockout only. The difference between the two results is originated by the 
imaginary part of the optical potential. 

FSI effects are similar on the longitudinal and transverse components of the 
nuclear response and are important both in the relativistic and in the nonrelativistic 
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calculations. The role and relevance of the various effects can however be different in 
the two frames. The final effect is similar and produces qualitatively similar results 
in comparison with data. The longitudinal response is usually well reproduced, 
while the transverse response is underestimated. This seems to indicate that more 
complicated effects, e.g., two-body meson exchange currents, have to be added to 
the present single particle approach. 
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Figure 1: Longitudinal (upper panel) and transverse (lower panel) response 
functions for the 12 C(e, e') reaction at q = 400 MeV/c calculated in the relativistic 

frame. Solid and dotted lines are obtained with and without the factor \Jl — V'(E), 
respectively. Dashed lines give the result without the integral in Eq. (31). Dot- 
dashed lines are the contribution of integrated single nucleon knockout only. The 
data are from Ref. [21]. (from Ref. [13]) 
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Figure 2: Longitudinal response function for the 12 C(e,e') reaction at q = 400 
MeV/c. Solid dotted and dashed lines arc the same as in the upper panel of Fig. 1. 
Dot-dashed and long-dashed lines are the nonrelativistic results with and without 
the factor Jl — V'(E), respectively. Data as in Fig.l. (from Ref. [13]) 
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Figure 3: The same as in Fig. 1, but for q = 500 MeV/c. The data are from Ref. 
[21]. (from Ref. [13]) 
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Figure 4: Longitudinal (upper panel) and transverse (lower panel) response 
functions for the Ca(e, e') reaction at q = 450 MeV/c. The Saclay data (open 
circles) are from Ref. [22], the MIT-Bates (black circles) are from Ref. [23]. Line 
convention as in Fig. 1. (from Ref. [13]) 
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